Observing the evaporation transition in vibro-fluidized granular matter 
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Abstract 

By shaking a sand box the grains on the top start to jump giving the picture of evaporating a sand bulk, and a gaseous 
transition starts at the surface granular matter (GM) bed. Moreover the mixture of the grains in the whole bed starts to move 
in a cooperative way which is far away from a Brownian description. In a previous work we have shown that the key element to 
describe the statistics of this behavior is the exclusion of volume principle, whereby the system obeys a Fermi configurational 
approach IJ. Even though the experiment involves an archetypal non-equilibrium system, we succeeded in defining a global 
temperature, /3~^, as the quantity associated to the Lagrange parameter in a maximum entropic statistical description. In fact 
in order to close our approach we had to generalize the equipartition theorem for dissipative systems. Therefore we postulated, 
found and measured a fundamental dissipative parameter, 5, written in terms of pumping and gravitational energies, linking 
the configurational entropy to the collective response for the expansion of the centre of mass (cm.) of the granular bed. Here 
we present a kinetic approach to describe the experimental velocity distribution function (VDF) of this non-Maxwellian gas of 
macroscopic Fermi-like particles (mFp). The evaporation transition occurs mainly by jumping balls governed by the excluded 
volume principle. Surprisingly in the whole range of low temperatures that we measured this description reveals a lattice-gas, 
leading to a packing factor, pf, which is independent of the external parameters. In addition we measure the mean free path, 
L, as a function of the driving frequency, /e, and corroborate our prediction from the present kinetic theory. 



The granular matter bed was set up with 5437 balls 
of Zr02 — Y2O3, in a glass container of Dc — 50mm 
diameter (which is open on the top) in a chamber at 
latm of air under low water vapor content conditions 
(5.8 ± 0.2gr/m^), see Fig. l(a,b). The diameter of the 
balls was D — 1.99mm and their mass m — 26.8±0.1nigr. 
Under vertical vibration, it is possible to get the balls 
at the lower gravitational position of the container in a 
crystalline state (the bulk), and those at the upper one 
in a gas state, which is called the fluidized gap |l|. This 
non-Markovian gas doesn't expand to the whole volume 
of the container, and the number of the active moving 
particles, in the gas phase, is a function of the pump- 
ing energy. Thus L of the mFp determines the occupied 
volume of the evaporated gas. A sinusoidal vibration is 
driven by a vibration plate on the GM bed of height 
h (with intensity F = AuP'jg, where A is the ampli- 
tude, g is the acceleration of gravity, and lo = 27r/e). 
The vibration apparatus is set up by an electromag- 
netic shaker (TIRAVIB 5212) which allows for feedback 
through a piezoelectric accelerometer for the control of 
/e and F in the range of 10-7000Hz, and 2 - 40.g re- 
spectively. The control loop is completed by an Oscilla- 
tor Lab- works SC121 and a TIRA 19/z amplifier of Ikw. 
The vertical trajectory (stochastic realization) z{t) of one 
particle was followed in a window of 12mm with a laser 
device by using a triangulation method, see Fig. 1(c) 
and its spectrum is shown in Fig. 1(d). A laser emitter 
with a spot of 70/im and a linear image sensor (CCD-like 
array) enables a high speed measurement with 100/isec 
sampling. The linear image sensing method measures the 
peak position values for the light spots and suppresses the 
perturbation of secondary peaks, which makes possible a 



resolution of l/j,m. The shaker and the laser displace- 
ment sensor were placed on vibration-isolated tables to 
decouple them from external perturbations. z{t) were 
taken with a 9354 C Le Croy Oscilloscope of 500MHz. 
The velocity, V{t) = dz/ dt, of the mFp was calculated 
numerically for At — 100/isec from z{t) registers. The 

amplitude dispersion ct^ = {z{t)^) — {z{t))^ and ve- 
locity variance cr^ = {V{tf) - {V{t)) were obtained 
from a window of 2 seconds for each pair of registers 
{z(i), ]/(<)}. In Fig. 2(a) we report cr^ against /e (from 
40-lOOHz) for fixed F = 10, for a GM bed of /i = 18mm. 
We test experimentally in Fig. 2(b) the relation between 
Gz and ay. Beyond the Gibbs ensembles, it should be 
possible to conceive a common conceptual framework un- 
der the non-equilibrium thermodynamics that involves 
jammed states, as described by Edwards ensembles 
From a granular gas it is possible to get either a jammed 
state, for example by reducing or quenching F at a given 
/e _,5i, or a crystal by slowly scanning up /e in a fixed 
accelerated experimental set-up (cooling down). In fact 
it is possible to see how a non-Brownian gas is condensed 
into a fee crystal 0] . 

The non-equilibrium GM temperature is a controver- 
sial quantity presenting anisotropy and inhomogeneous 
characteristics, and its behavior as a function of the ex- 
ternal pumping of energy is yet an unsolved question. 
However we have laid out an experimental set-up that 
makes it possible to measure the "global" temperature 
as the Lagrange parameter that appears in 0]. 

Since energy is constantly being added to the system 
a nonequilibrium steady state (s.s.) mass profile for the 
GM bed is reached. We report here, upon measuring 
the evaporation transition line in a weakly vibrated GM 




FIG. 1: (a) Top view of the GM bed in a condensed fee state, 
(b) Set up for distance measurement (amplitude realizations 
z(t)) corresponding to 11 layer GM {h — 18mm) with Zr02 — 
Y2O3 balls {D — 1.99mm) in a glass container of Dc — 50mm 
diameter, (c) Stochastic realization z{t) at 83Hz, lOg, and (d) 
the corresponding anomalous power law spectrum with 
v = 3. 




FIG. 2: (a). Evaporation transition for the experiment of 
Fig.l. az as a function of the external driving frequency f^. 
(b) az vs. ay and the corresponding theoretical fittings. The 
fe was from 50 — lOOHz. 



medium, the dramatic consequences of the exclusion of 
volume principle, for a typical non-equilibrium GM sys- 
tem. It is important to remark that below this criti- 
cal point the GM behaves as a crystal where the parti- 
cles move in "thermally" activated harmonic oscillations, 
thus in this regime a Brownian description is suitable, 
and this fact was already reported by 0|- In a previ- 
ous work we have studied the spectrum of vibrated 
GM under gravity, and shown that in the weakly excited 
regime the dynamics of these evaporated particles pre- 
cludes describing them as Brownian. They show (over 
2 decades) a power law spectrum Sz{f) ~ 1//" with 
an anomalous exponent v {2 < < 5) that depends on 
the excited collective movement of the grains. This fact 
leads us to the conclusion that a description of the present 
chaotic cooperative dissipative dynamics of the GM par- 
ticles, must be done in terms of non-Markov particles 
representing the collective effects of the excited GM0,0. 
Recent reports concerning the non-Brownian behavior of 
the velocity fluctuations can also be found in resistivity 
experiments which show that the electric noise has inter- 
esting properties of scale invariance and of intermittency 
which arise from thermal expansion locally creating or 
destroying electrical contacts addition, associ- 

ated anomalous transport properties have been reported 
as non-Brownian motion jlll] . 

Following Hayakawa and Hong |3| we have introduced 
an approach to understand the behavior of a weakly ex- 
cited GM Jj]. Our 3— dimension A^— particles experimen- 
tal conditions allow us to consider a system of rows 
in a cylindrical container as a 1 dimensional degenerate 
Fcrmi-likc system. Our laser facility was set up to inves- 
tigate the occupation dynamics at the fluidized gap that 
appears at the top of the n-rows GM bed. In Fig. 1(c) a 
typical z(t) is shown. This measurement corresponds to 
the beginning of evaporation, where the spectrum expo- 
nent is u = 3, thus we shown that the velocity fluctuation 
does not correspond to the one that could be obtained 
from a Brownian particle 6] . 

By focusing on the configurational properties of an ex- 
cluded volume theory the s.s. mass profile can be un- 
derstood in terms of a configurational maximum prin- 
ciple assumption. Excluded volume interactions of the 
GM do not allow two grains to occupy the same state 
of gravitational energy. Following Landau's approach 
to study a non-equilibrium system, the maximization of 

= InJli [^^-/N.lin - NiY.] yields that the density pro- 
file is 4>{e) = [1 + Q exp{l3e)]~^ il/ mgD where il is the 
degeneracy (number of boxes in which to put the maxi- 
mum amount of balls in a constant gravitational layer); 
/3 is a Lagrange multiplier parameter, Q^^ = cxp(/3/i) 
the fugacity, and e — mgDs with s = 0, 1, 2, 3 • • • . In- 
troducing the normalization we get the relation between 
the chemical potential /j, and the temperature , i.e., 
eyij>{NPmgD/D,) = l + exp(/3/i). The zero-point "chemi- 
cal potential" is /io — mgDN/ft, so N/il is the number of 
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balls in an elementary column of diameter D. It is trivial 
to understand that without pumping of energy the cm. 
is characterized by Zc.m. = fJ-o/'^mg = h/2. For a weakly 
vibrated GM we have shown that /3 is a non-trivial func- 
tion of the ay 0] (if the particles were Brownons the 
relation should be | v |(x /3^^/^). 

Shaking gently a box of GM is an example of a nonequi- 
librium system that can be characterized by a global dis- 
sipative parameter S, which in fact is given as the ratio 
of the two important energies of the system, 



mgAzc 



Alo 



(1) 



Since the motion of these mFp can be studied in terms 
of the cumulative probability — mgz) / N, z = Ds 
we have shown that at low temperature, /3/io S> 1, we 
get Q 



a, ~ 2e/mgl3. 



(2) 



where 6 = ln(l + y^l - Ij^fe) + \. This expression 

gives the amplitude dispersion = \J (z{t)^) — {z{t))'^ ~ 
a^/mg as a function of [3. The connection between the 
fluctuation in the velocity cr^ = {V{tf) - {V{t))'^ and (3 
is given by our generalized equipartition law: ^m*ay — 
i^^/3~^, where m* /m = S and AN/ N is a relative fac- 
tor that counts the "active" mFp in the fluidized gap. 
From this we get 



l(^ln[2e^--l]-l)=W^, 



(3) 



and this formula connects ay, the dissipative parameter 
5 and f3. So using ^ we relate the collective response of 
the system: the expansion of the cm. Azc.m. = U/N — 
/^o/2, where U — Jq ^'t'i^) de, with the variance ay and 

fe- ^ 

Our experiments show that the fluidized gap reveals 
an inelastic gas heated in a non-uniform way. Now let 
us understand this gas of mFp by presenting an effective 
kinetic theory. The collision frequency between mFp is 
given by = Vagn, where cto = ttD^ is the total scat- 
tering cross section between two hard spheres, V their 
mean relative speed, and n is the mean number of these 
mFp per unit volume. The mean free path is L = tv. 
Proposing the relative speed V = ciay and the mean 
speed V = c^ay (actually V should be somewhat larger 
than v) and using that the mean number of "active" mFp 
per unit of volume is given by 
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FIG. 3: (a). Mean-free-path L, and (b) the Lagrange param- 
eter /3 vs. frequency /e. 



and assuming all the constants ci = C>{1) we get for 
/3/xo » 1 



2\Dl AN' 



(5) 



L calculated from the experimental data {a^, AN (ay)} 
is depicted in Fig. 3(a) as a function of fe, in addition in 
Fig. 3(b) we show /3 against fe to note the evaporation 
transition point. From the linear fitting of Fig. 3(a) we 
obtain L ~ 2.28mm, this value is indeed only an esti- 
mation since we don't know the values of the constants 

With a video camera microscope we took photographs 
of the surface of the GM-bed from the top of the con- 
tainer and we found that time average pictures show that 
L ^ D, this result agrees with the measurements of cor- 
relation in the gas phase estimating the degen- 
eration fl in the expression Q from Tr{Dc/2) ^ flD^, 
we get L ~ AQD/TT\n2 ~ 1.351?. This important result 
shows that a "jumping model" can therefore be consid- 
ered to build up a kinetic transport theory; thus allowing 
to calculate the VDF for the mFp in the fluidized gap. 
We want to emphasize that Fig. 3(a) supports the ap- 
plication of a jumping model to study the kinetics of this 
non-Maxwellian gas because L ^ D, i.e., we shall asso- 
ciate a jump probability distribution P{r) (proportional 
to ipipo + e) l^f I J where r = e/mg > 0, [l|) with the 
transport mechanism of a set of independent particles in 
a background of a well defined global temperature (3~^. 
A crucial ingredient for understanding the dynamics of 
the GM would be to be able to build up a statistical me- 
chanics for the construction of the VDF of the driven GM 
system. A less ambitious program is to have a model for 
the kinetic transport of the mFp and that is what we are 
going to present. 

Supposing that during the time between collisions the 
velocity V of each particle is essentially constant, a jump 
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of length r is performed with velocity V = r/r > 0; then 
defining the kinetic energy as K = mV'^/2 and using the 
transformation of random variables we get P{K)dK — 
P{r)dr. This "jumping model" allows us to arrive at a 
closed expression for the density of the kinetic energy in 
terms of the /3 of the s.s. of the configurational entropy, 
i.e., 

P{K)^gTJ^i,{fio+gTV2^); / P{K) dK ^ 1. 
V 2J< Jo 

(6) 

Using the mass profile we get ipi^e) = —d(j)/de, there- 
fore from (jSJ wc obtain an analytic formula for the ki- 
netic distribution function (KDF) P{K), which predicts 
at medium and high energy ranges a larger population 
than the Maxwellian. Note that in the asymptotic behav- 
ior e » ^0, (ie~^'^ then KDF for high energy goes 
like ~ gTyJ m/K g-a'^P^/^rnK jg non-Maxwellian. 

An important program is to perform a suitable exper- 
iment to be able to determine the functional form of the 
KDF, and this is what we are going to present in the 
rest of the paper. While in the first experiment we scan 
in fe at a given F to investigate the behavior of /3, in a 
second experiment we carry out an isothermal-like mea- 
surement during a long enough time to get the KDF and 
the corresponding /?. 

From the measurement of z(t) we were able to calcu- 
late the time dependence of the kinetic energy fiux in 
the measurement of a 1-dimensional window, from which 
it is possible to obtain the KDF. At a given /e and F 
we measure 5000 windows of 2 seconds to get a total of 
25 millions of events from which we build up the s.s. 
KDF. In Fig. 4 we fit our experimental data by us- 
ing the formulae ©; the experiment was carried out for 
fixed fe — 80Hz and F = 10. Medium and high energy 
particles, which are at the surface of the vibrated GM 
bed, are well described by formulae while at low en- 
ergy the behavior is Maxwellian. Also we fit the whole 
experimental results with a Maxwellian distribution in 
kinetic energy ll{K) — exp {—Kf3) / yJ-nKj fi. As well as 
with a stretched exponential for the velocity distribution: 
n(y) ~ exp [- (V"/Vbrj, with a = 3/2, where Vb is the 
thermal rms velocity [l^; this last fitting shows that the 
agreement is also good. Since Olafsen et al. reported 
the stretched exponential law for the VDF of the gas 
phase of a vibrated monolayer GM^^, several research 
gro ups have been measuring this law for granular gases 
|l5j . A delicate obstacle for the analytical solution of 
this problem is the supply of energy to the granular bed 
which compensates the dissipation caused by inelastic col- 
lisions [l^. an aspect which was clarified by van Zon el 
at. using simulations of molecular dynamics For 
granular gases there is indeed a family of distributions, 
with apparent exponents a < 2 which are governed by 
the ratio between the numbers of heating events and in- 
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FIG. 4: Experimental kinetic energy distribution function 
(KDF) and tlie corresponding fittings according to our "jump- 
ing" kinetic model, the Maxwellian and the stretched expo- 
nential cases. 

elastic collisions |l6| . Either under uniform heating or 
boundary heating of the GM, if the dissipative collisions 
dominate the heating process produces a liquid-like clus- 
ter surrounded by a gas phase and the VDF is strongly 
non-Gaussian. Also after Olafsen et al. the conclusions of 
several authors keep clear the conceptual relationship be- 
tween this high energy tail behavior, the exclusion of vol- 
ume and the clustering in granular gases. Our description 
adds a new picture for clarifying this delicate point, the S 
parameter defined in equation JQ) and the generalization 
of the equipartition law Q permit us to tackle this obsta- 
cle, to give the KDF in terms of /3. Our approach was ob- 
tained from the physical parameters: {fe, F, /ip, D, N, ft} 
and by using the free parameter S from the experiment 
described in Fig. 2 and 3. The fitting for the curve in Fig. 
4 gives f3 = (1.7 ± 0.5)(/iJ)~^ and for the corresponding 
value from Fig. 3(b) we get (5±l)(/xJ)^^, which in fact is 
a good agreement considering that in the first experiment 
we are scanning in fe at constant rate, and in the sec- 
ond experiment (time average) both fe and F are fixed. 
The high energy tail corresponds to the beads near to 
the surface of the GM bed where the jumping combined 
with collisions dominates the process. 

A point of special interest is that the packing fraction 
for this fluidized gap is constant for the whole range of 
low temperatures we used, which can be estimated as: 

2 / dV D-AN 7rln2 

At constant cooling down by increasing fe at constant 
F, this lattice-gas condenses into an fee crystal, and vice- 
versa from the fee lattice under constant heating rate the 
expansion of the cm. entails the expansion of this fee. 
This expansion under gravity creates fluctuations in the 
number of particles n(e) around e = /io, which induces 
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the transport of matter at the fluidized-gap. The time 
integration of such fluctuations (as we did during 3hs. in 
0) reveals that the fluidized-gap has a structure with 
a Fermi profile. Thus we conclude that for our experi- 
mental conditions the gaseous gap is indeed a lattice-gas 
{fee where half of sites are occupied) with pf according 
to equation ((T)). Nevertheless, the picture we get for this 
lattice-gas should be understood as a "jammed state of 
motion", since the pf does not change for the range of 
low temperatures we measured. 

The important experiment of D'Anna et al. shows 
that a vibrated granular bed exhibits a formal anal- 
ogy with a thermal bath, where it is possible to apply 
a fluctuation-dissipation theorem. In that experiment a 
harmonic oscillator of frequency too (the torsion pendu- 
lum) is coupled to the Brownian particles in the vibrated 
granular bed. This thermal bath is generated by driv- 
ing the GM in a frequency band between 300 — 900Hz at 
a given F (from 1 to 11. 651); in fact, out of equilibrium 
this thermal bath is thought within the frame in 21 . The 
D 'Anna's torsion pendulum is this thermometer coupled 
to the Brownian ensemble, the driven fluctuation pumps 
energy while its interaction with the same Brownons pro- 
vokes the dissipation of its energy. Nevertheless if evapo- 
ration starts the viscosity of the thermal- like bath should 
therefore be related to the exclusion of volume, since the 
jumping of the particles is conditioned by the availability 
of free volume. 

In a system in equilibrium and in D 'Anna's exper- 
iment, equipartition means that K must be kBT/2 
and by measuring long enough they get a temperature 
for the thermal bath (or granular bed), the range of 
temperatures tested were lower than O.l^J. Our laser- 
thermometer explores beyond the Markovian range of the 
GM bed, which happens when the jumping events start 
(evaporation) and the fluidized gap appears. Our ther- 
mometer is a test particle of the GM bed which occupies 
the position tested by the laser. 

Measurement of /?"^ as a function of /e exhibits at 10 
PKelvin the critical point that characterizes the evapo- 
ration transition. For /3~^ < lOPKelvin (/3~^ < O.l^J) 
the system is condensed and the motion corresponds to 
Brownian oscillators. If evaporation starts the following 
picture is possible: when our ensemble of mFp is per- 
turbed, the energy is dissipated by the mechanisms of 
jumping (under gravity) and random inelastic collisions 
between the mFp, but conditioned to the fact that these 
events are self-organized in time, keeping constant the 
Pf oi the fluidized gap. The relation between the trans- 
port and the number of "active" particles AN (at the 
Fermi sea) links the two Lagrange parameters: (3 and 
H (compactivity) , and therefore describes the existence 
of a lattice-gas. The coupling between /9~^ and fj, tells 
us that the important rule, for the ensemble of mFp, is 
to dissipate energy keeping constant the lattice struc- 
ture. We remark that the difference with the well known 



Markovian motion is that Brownons can dissipate energy 
through viscosity independently of any gas structure. 

Our approach focuses on the potential energy which 
is a function of the configurational entropy and the re- 
lation with the kinetic energy per active particle at the 
gap. In other words, while in the Brownian situation 
we can describe the system in terms of energy fluctua- 
tions, viscosity and temperature, beyond the Markovian 
approximation, the important events are the fluctuations 
in the number of particles. Then we are indeed report- 
ing on the existence of a eompensation relation between 
the fluctuations in the number of mFp at the gap (or 
the matter transport), the Fermi profile held by the S 
dissipative parameter, and the global temperature /3~^. 
A fluctuation in the number of mFp in the gap will be 
compensated in time in such a way that the lattice-gas 
keeps Pf constant: this is a fluctuation-compensation re- 
lation. Therefore a stochastic formalism is still valid but 
the equations should describe a fluid where the complex- 
ity is the description of the lattice-gas, the jumping phe- 
nomena and collisions in a space framed by the exclusion 
of volume. 

This approach agrees with the discussion in |3| related 
to the idea that systems far from equilibrium can be de- 
scribed with more than one temperature. Such tempera- 
tures are related to time constants of the different config- 
urational structures. For a glass, the time for relaxation 
to equilibrium is extremely long, and it is possible to de- 
scribe this system by using more than one thermometer 
with different time constants. In GM experiments, the 
structure of the gap is only measurable if we wait enough 
time. This structure is related to transport of mass at 
the gap, as is revealed by our generalized equipartition 
law. The characteristic time we must wait to measure /3 
in this system is the time we need to integrate the mass 
profile. However unlike glasses, GM vibrated under grav- 
ity is far away of equilibrium but stationary: the system 
exists while the S parameter is held at a fixed value, and 
the fluctuation-compensation should play a role stabiliz- 
ing the steady state, i.e., the precise balanee of fluctua- 
tions permits us to know the global temperature without 
waiting too long. 
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